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ABSTRACT 

We investigate the production of electron-positron pairs by inverse Compton 
scattered (ICS) photons above a pulsar polar cap (PC) and calculate surface heating 
by returning positrons. This paper is a continuation of our self-consistent treatment 
of acceleration, pair dynamics and electric field screening above pulsar PCs. We 
calculate the altitude of the inverse Compton pair formation fronts, the flux of 
returning positrons and present the heating efficiencies and X-ray luminosities. We 
revise pulsar death lines implying cessation of pair formation, and present them in 
surface magnetic field-period space. We find that virtually all known radio pulsars 
are capable of producing pairs by resonant and non-resonant ICS photons radiated 
by particles accelerated above the PC in a pure star-centered dipole field, so that 
our ICS pair death line coincides with empirical radio pulsar death. Our calculations 
show that ICS pairs are able to screen the accelerating electric field only for high 
PC surface temperatures and magnetic fields. We argue that such screening at ICS 
pair fronts occurs locally, slowing but not turning off acceleration of particles until 
screening can occur at a curvature radiation (CR) pair front at higher altitude. In the 
case where no screening occurs above the PC surface, we anticipate that the pulsar 
7-ray luminosity will be a substantial fraction of its spin-down luminosity. The X-ray 
luminosity resulting from PC heating by ICS pair fronts is significantly lower than the 
PC heating luminosity from CR pair fronts, which dominates for most pulsars. PC 
heating from ICS pair fronts is highest in millisecond pulsars, which cannot produce 
CR pairs, and may account for observed thermal X-ray components in the spectra of 
these old pulsars. 
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1. INTRODUCTION 

In the last several years, the basic model of particle acceleration above a pnlsar polar cap 
(PC) has been undergoing significant revision. Sturrock (1971), Ruderman & Sutherland (1975) 
and Arons & Scharlemann (1979) originally proposed that particles are accelerated by an induced 
electric field, producing curvature radiation (CR) photons which create electron-positron pairs 
in the strong magnetic field. The pairs short out the electric field above a pair-formation front 
(PFF), self-limiting the acceleration. In the process of screening or shorting-out the electric field, 
some fraction of the positrons is accelerated back toward the PCs and heats the surface of the 
neutron star (NS), producing a potentially observable X-ray emission component (Arons 1981). 
Two recent developments have introduced important changes to this picture. First was the finding 
of Muslimov & Tsygan (1992) that the effect of inertial-frame dragging near the NS surface 
greatly increases the induced electric field E \| above the PC in space-charge limited flow models 
(Arons & Scharlemann 1979). The second was the realization that inverse-Compton radiation of 
the primary particles can produce pairs which could potentially screen the electric field (Zhang et 
al. 1997, Harding & Muslimov 1998; hereafter HM98). 

We have investigated the electric field screening and PC heating in the revised space-charge 
limited flow (SCLF) model. In the first paper (Harding & Muslimov 2001; Paper I), we have 
presented our results for screening and PC heating by CR PFFs, as in Arons (1981), but using 
the frame-dragging electric field of Muslimov & Tsygan (1992). Paper I outlined a self-consistent 
calculation of the PFF height, returning positron flux and the screening scale length, where the 
E || and the primary flux are adjusted for the change in charge density caused by the returning 
positrons. A main assumption of this calculation was that the screening scale is small compared 
to the height of the PFF, the location where the first pairs are produced. This turns out to be 
generally true for screening by CR-produced pairs because the pair cascade multiplicity grows 
rapidly over small distances due to the strong dependence of CR photon energy on particle energy. 
Thus, the existence of a CR PFF always results in PC heating and full screening of E\\. We 
found that the PC heating luminosity, as a fraction of the spin-down luminosity, increases with 
pulsar characteristic age, r = P/2P, and should be detectable for pulsars with r 10 6 yr. The 
most significant heating occurs for pulsars near the death line for CR pair production, which is 
t ~ 10' yr for normal period pulsars and r ~ 10 8 yr for millisecond pulsars. Our predicted X-ray 
luminosity due to PC heating is about a factor of ten higher than the X-ray luminosity predicted 
by Arons (1981), due to the increase in accelerating voltage drop resulting from the inclusion of 
inertial-frame dragging effects. We also predicted that older pulsars should have higher PC surface 
temperatures from heating. 

In this second paper (Paper II), we present results of our investigation of electric field 
screening and PC heating by inverse-Compton scattering (ICS) radiation PFFs. We investigate 
PC heating by ICS produced pairs in all pulsars, including those which produce CR pairs. For the 
older pulsars that do not produce pairs through CR, PC heating by positrons from ICS cascades 
is especially important. Our treatment of the ICS PFFs follows closely that of HM98, where 
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cyclotron-resonant and non-resonant scattering are considered as separate radiation components. 
HM98 found that the ICS PFFs are located much closer to the NS surface than are the CR PFFs, 
because primaries with Lorentz factor of only 10 4 - 10 5 can produce pairs via ICS whereas 
Lorentz factors of at least ~ 10 7 are required for production of pairs via CR. They also found that 
the ICS PFFs of positrons returning from the upper PFF occur at a significant distance above the 
surface, so it is possible that ICS pairs fronts are unstable if screening of E y at the lower PFF 
occurs. However, HM98 did not determine whether the ICS pairs were capable of screening either 
at the upper or lower PFFs. 

This paper will attempt to answer those questions with a calculation similar to that of Paper 
I. In Section II, we describe the acceleration model for this calculation. We use a solution for E\\ 
from Poisson’s equation which imposes an upper boundary condition requiring that E^ = 0 at 
the location of screening (as in HM98), which differs from that of Paper I which used a solution 
with no upper boundary. The solution with an upper boundary is more appropriate for screening 
at ICS PFFs which form close to the surface, where the screening scale is generally comparable 
to the PFF height. This solution is required here because pairs from ICS photons may screen 
E || close to the NS surface. In this case the value of E y is suppressed due to proximity of the 
screening to the surface. Section II also presents a calculation of ICS PFF height and the location 
of a new ICS pair death line. We find that virtually all known pulsars can produce pairs by 
either resonant ICS (RICS) or non-resonant ICS (NRICS) in a pure dipole field. Hibschmann & 
Arons (2001) (hereafter HA01) have recently taken a different approach to determining pulsar 
death-lines including ICS produced pairs. They define the PFF as the location where the pair 
multiplicity achieves that required for complete screening of the En, whereas we define the PFF as 
the location where pair production begins (i.e. where the first pairs are produced). In the case of 
CR, this distinction is minor since the screening scale is small compared to the PFF height, but 
in the case of ICS the distinction is very important since the screening scale is comparable to the 
PFF height. HA01 also assume that the pair multiplicity required for screening is determined by 
the difference between the actual charge and the Goldreich-Julian charge at a distance from the 
NS surface that is roughly equal to the PC radius, whereas in our calculation the charge density 
required for screening is determined by the charge deficit at the location of creation of pairs which 
can be much smaller. With their more restrictive definition of the PFF, their death lines differ 
significantly from ours. In Section III, we give self-consistent solutions for the fraction of positrons 
returning to the PC and for the screening scale. We also explore the question of whether positrons 
returning from the upper PFF can screen E\\ near the lower PFF close to the NS surface. The 
subject of pulsar pair death lines, i.e. which NSs are capable of pair production, will be discussed 
in Section IV. In Section V, we present our calculations of PC heating luminosity due to returning 
positrons, giving both numerical and analytic estimates. Summary and conclusions, as well as a 
comparison of results from both Papers I and II will be given in Section VI. 
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2. THE ACCELERATION OF PRIMARIES AND ONSET OF PAIR 

FORMATION 

2.1. Acceleration model 

In this paper we exploit the same model for charged particle acceleration that is described in 
Paper I. Namely, in the acceleration region we use the appropriate solution for the electric field 
and potential presented in Paper I. For example, for most ICS pair fronts, where the screening 
occurs over the length scale smaller than the PC size, we use our formulae (A7)-(A10) of Paper I, 
whereas in the screening region (see § 3.2 for details) for all cases we model the electric field by 
equation (26) of Paper I. 


2.2. Altitudes of the pair-formation fronts 

2.2.1. Analytic estimates 


For our analytic estimates we shall use the same simplified expressions for the accelerating 
electric field (cf. eqs. [34], [35] of Paper I), 


17 Rl2 

E \\6 = - p “ K 0.15 


1.3 z P~ 1/2 
5 • id -3 p- 1 


(i) 


where the upper expression corresponds to the unsaturated regime (rising part of the accelerating 
field), and the lower - to the saturated regime (nearly constant accelerating field). Equation ([|) 
assumes that £ = 0.5 and cos x ~ 1; Ey 6 = Em/ 10 6 esu, K0.15 = «/0.15, B 12 = -Bo/10 12 G, and 
z is the altitude in units of stellar radius, where £ is the magnetic colatitude in units of the PC 
half-angle, x is the angle between the magnetic and spin axes, Em is the component of the electric 
field parallel to the magnetic field, Bq is the surface value of the magnetic field strength, k is 
the dimensionless general relativistic parameter originating from the frame-dragging effect and 
accounting for the stellar compactness and moment of inertia. 

Note that throughout this paper in all practical formulae we assume that P is dimensionless 
value of pulsar spin period measured in seconds. Also, in all our analytic estimates based on the 
above formula we shall discriminate between unsaturated and saturated regimes of acceleration, 
with the unsaturated/saturated regime occuring in the case where the characteristic altitude of 
pair formation is smaller/larger than the PC radius. The formal criteria corresponding to these 
regimes, e.g. in B-P diagram, depend on the radiation mechanism for pair producing photons 
and can be derived from the conditions £* < 1 and £* ~ 1 for unsaturated and saturated regimes, 
respectively, where £* is the characteristic altitude of screening scaled by the PC radius (see eqs. 
[|36|]-[|38|] below). These criteria translate into 


P < P* Unsaturated, regime 


(2) 
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and 


where P* is defined as 


P r^i P * Saturated regime , 


( 3 ) 


Curvature radiation 

P* (Ci?) = 0.1 Bti 9 - (4) 

Before we present P* for the ICS case, note that for the physical conditions we discuss in this paper 
it is convenient (especially in our analytic calculations) to treat the ICS photons as generated 
via two different regimes of scattering: resonant (R) and non-resonant extreme Klein-Nishina 
scattering (NR). In reality, the spectrum of ICS photons from an electron of Lorentz factor 7 , a 
subset of which produce pairs, may be produced by scatterings in both the R and NR regimes. 
However, the first pairs which mark the PFF location will come from one of the two regimes. 
Thus, in our analytic calculations, to avoid unnecessary complexity, we shall differentiate between 
R and NR regimes. Now let us present the expressions for P* corresponding to these regimes 

Resonant ICS 

P* (i?) = 0.1 B% 7 , (5) 

Non-resonant ICS 

Pi NR) = 0.4 P^ 7 . ( 6 ) 


In this Section we shall explore the altitudes of PFFs produced by the ICS photons. To 
estimate the altitude of the PFF above the stellar surface we can use, as in our previous papers 
(see HM98, eq. [1]; and Paper I, eq. [37]), the following expression 

So = min[S' a ( 7 min ) + *Sp(^min )], (7) 

where >S , a (7 in in) is the acceleration length that is required for an electron to produce a photon of 
energy £ m i n , and S p (s min ) is the photon pair-attenuation length. 

Resonant ICS 

In the R regime of ICS the characteristic energy of a scattered photon is 

£ ~ 2 7 P', (8) 

where B' = B/B cr (< 1) is the local value of the magnetic field strength in units of critical field 
strength, B cr = 4.41 x 10 13 G. 

In this paper in our analytic estimates we assume that B ~ Bo, which is justified for all cases 
except for the case of millisecond pulsars. By substituting s into the expression for S p (see HM98 
or Paper I for details) we get 

S (R) = cf V 7 , (9) 

where c\ R ^ = 2.2 x 10 10 P 1 // 2 /P ^ 2 cm ■ 
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As far as the acceleration length S a is concerned, it is independent of the radiation mechanism 
(provided that the radiation losses are negligible during acceleration), and the formulae of Paper I 
(see eqs. [38]) are still applicable. After substituting expressions for S a and Sp ^ into © we find 
that Sq is minimized at 


7 


(R) 


= 10 b 


0.9 P- 1 / 6 ^ 1 if P<P* (R) , 
0.2 p- 3/4 s- 1/2 if p A p( R) . 


( 10 ) 


Now we can evaluate expression © at 7min to calculate the dimensionless altitude (scaled by 
a NS radius) of the PFF 


JR) = Q (R) f JR) 


= S^W/R = 10 


-2 


7 P 2 / 3 B{ 2 1 
17 P 5 / 4 P“ 3/2 


(i?) 


if P < P| 
if P ~ Pi 


( 11 ) 


Non-resonant ICS 

In the NR regime we can use the extreme Klein-Nishina formula, and write for the energy of 
the scattered photon 


e ~ 7 . 


( 12 ) 


The corresponding photon pair-attenuation length is 

S (NR) = c (NR) h ^ 

where c\ NR ^ = 10 9 P 1//2 /Pi2 cm. The expressions for 7^ r f^ and Zq NR ^ then read 

if P<P* (7VR) , 
if P~ P^ NR \ 


cm = 10 s 

/ min 


p- 1 / 6 p 1 " 2 1/3 

0.6 P“ 3 / 4 


(13) 


(14) 


and 


SNR) 


= 10 


-2 


3 (P/Pi 2 ) 2/3 if'P<P (jVR) 

4 P 5 / 4 /P 12 


if p~p! nr) . 


(15) 


The unsaturated regime of formula (14) agrees within a factor of 
calculations for the ms and middle-aged pulsars. 


2 with our numerical 


2.2.2. Numerical calculations 

The location of the PFF for the different pair-producing radiation processes is more accurately 
determined using numerical calculations of the minimum height at which the first pairs are 
produced. Here we incorporate full spatial dependence (in r and 9) of quantities such as the 
magnetic field and radius of curvature in computing the altitude of the PFF. However, as we will 
discuss in Section 3.2, we use a one-dimensional model for the pair dynamics to treat the screening 
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of the electric field about the PFF. The numerical calculation of the altitude of the PFF follows 
closely that of HM98, but with one notable exception. As discussed above, we use the solutions for 
Eu with an upper boundary to compute the ICS PFFs (as in HM98) for normal pulsars, whereas 
we use the solutions for Eu with no upper boundary (but which saturates at infinity) to compute 
the CR PFFs (as in Paper I) and the ICS PFFs for ms pulsars. Otherwise, we have used the same 
expressions for the energy of pair producing photons: equation (29) of HM98 for CR photons and 
equation (43) of HM98 for ICS photons. We have also made the assumption (as in HM98) that the 
primary electrons travel along the magnetic axis in computing the ICS energy loss and scattering 
rate. Unlike in the analytic estimates presented, we do not separate the RICS from the NRICS 
photons but take a weighted average for the ICS photon energy: 


{e (R)^(R) + {e W)rfNR) 

{ ' j(R) + 7 {NR) 


(16) 


Calculations of zq as a function of the colatitude £ for four different sets of pulsar parameters are 
shown in Figure 5. Comparison of Figures 5a and 5b confirms that the value of zq is independent 
of PC temperature, since RICS photons have average energies of 2y B' and NRICS occurs primarily 
in the extreme Klein-Nishina limit, where the average scattered photons energy is 7 . 


With numerical calculation of the PFF altitude, we are able to determine the pulsar 
parameter space in which the formation of a PFF by the different photon processes is possible 
(i.e. where 5p(e m i n ) < 00 ). Figure 1 shows the CR and ICS PFF parameter space as a 
function of pulsar period P and surface field strength Bq, as determined by the dipole formula 
Bq = 6.4 x 10 19 G (PP) 1 / 2 . We have also plotted observed pulsars from the ATNF catalog 

, which includes pulsars from the Parkes Multibeam Survey 
(Manchester et al. 2001) having measured P. For the Bq and P values above the lines in Fig. 

1 for the different processes, we were able to find numerical solutions to equation (|t|). Pulsars 
below the lines are therefore not capable of producing pairs from that process. We find that the 
majority of pulsars cannot produce pairs through CR. This result has been well known for some 
time (Arons & Scharlemann 1979, Arons 1998, Zhang et al. 2000). We find that virtually all 
known pulsars, with the exception of only a few ms pulsars, are capable of producing pairs via 
ICS. For the majority of pulsars, those with Ro~0.1R cr , the pairs are produced from NRICS 
photons. We have plotted ICS PFF boundaries for three different temperatures and, as expected, 
there is only a small dependence on PC temperature. The slight variation due to PC temperature 
occurs for higher field pulsars where scattering is not in the extreme Klein-Nishina regime and 
the photon energies become somewhat temperature dependent. The failure of the lowest-held 
ms pulsars to produce ICS pairs in our calculation may be due to the failure of our E y solution, 
obtained in the small-angle approximation, to accurately model the accelerating held in the ms 
pulsars which have large PCs with 9 pc ~ 0.3; or to our use of a canonical NS model in the case of 
ms pulsars which may have undergone signihcant mass accretion in Low-Mass Binary systems. 
We are currently exploring these effects to correct this shortcoming. Also, we shall examine the 
possibility of photon-photon pair formation in ms pulsars. We will discuss the broader implication 
of pair death lines in Section 4. 


(http://www.atnf.csiro.au~pulsar/ 




3. Returning positron fraction and electric field screening by pairs 


3.1. Analytic estimates 


To estimate the maximum fractional density of positrons returning from the upper PFFs 
we can use equation (33) of Paper I and the expressions for zq derived in § 2.2 . We must note, 
however, that in Paper I (where we investigated the CR case alone) in our analytic calculations 
of the fraction of returning positrons, we justifiably assumed that the screening occurs within 
the length scale much less than zq. This is not the case for ICS, for which the screening scale is 
determined by the scale of growth of energy of pair-producing photons. For the ICS photons the 
energy of pair-producing photons ~ 7 , so that the scale of growth of 7 governs the screening scale 
A s , i.e. 


As ~ 


7 


d'y/dz 


z=z 0 


0.5 Zq if P<P^ R ’ NR \ 

Z 0 if p>p( R ’ NR \ 


(17) 


Thus, in the case of ICS the screening effectively occurs over a more extended region. This means 
that even within the screening region the primary electrons keep accelerating and generating 
pair-producing photons. Since the pair-production rate per primary electron in the case of the 
ICS photons is not as high as in the case of CR photons, the assumption that the altitude of the 
screening is Zq would result in a significant underestimation of the returning positron fraction 
at the ICS-controlled pair-fronts. To come up with better estimates we should add A s to the 
calculated values of zq. or, according to formula (|l7|), multiply zq by a factor 1.5 and 2, for the 
unsaturated and saturated regimes, respectively. 


Evaluating expressions for p + /pcj (here is the charge density of returning positrons, and 
Pgj is the Goldreich-Julian charge density; see also eq. [33] of Paper I) at z* = zq + A s , we get 

Resonant ICS 


and 

Non-resonant ICS 



2P 2 / 3 /B l2 if P<Pi R \ 
hP 5 ^/B% 2 if 


(—) 

\PGjJ 


(NR) 


= 10 


-3 


7(P/B 12 ) 2 / 3 if P<P!) NR \ 

9 P 5 / 4 /R 12 if P ~ P* NRS> . 


(18) 


(19) 


In their paper HA01 assume that the fraction of returning positrons (see § 5.3) is of order of 
Kr pc /R, where r pc {K, 10 4 /v / P cm) is the PC radius, and R is the NS radius (= 10 6 cm). In 
contrast, in our calculations the fractional density of returning positrons is proportional to the 
dimensionless screening altitude (see Paper I), and can be roughly estimated as 1.5[re/(1 — k)]zo- 
Thus, the ratio of our value for the density of returning positrons to that of HA01 is of order of 
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So/r pc , which for the observed pulsar parameters may range from ~ 0.1 to ~ 20 . More specifically, 
for short-period pulsars our calculated fractions of returning positrons are up to a factor of ~ 10 
less than those of HA01, while for relatively long-period pulsars they are larger by about the same 
factor. Besides this quantitative difference, there is a significant difference in the dependence of 
our fractional densities on B and P. Note that in reality, the screening scale and thus also p + /pcj 
is dependent on the PC temperature, as will be shown in the numerical calculations. We cannot 
model this dependence analytically since we have assumed that the A s is a constant multiple of zq 
for all pulsars. 


3.2. Numerical calculations 


Our numerical calculation of the returning positron fraction and scale length of the ICS pair 
screening follows that described in Paper I for CR pair screening. As detailed in that paper, we 
first compute the pair source function in energy and altitude above the PFF and then compute 
the dynamical response of the pairs to the E\\ above the PFF by solving the continuity and energy 
equations of the pairs to obtain the charge density. The E\\ above the PFF is parameterized as 
an exponential with scale height A s , which is readjusted at each iteration to equal a multiple of 
the height at which the computed charge density from the pairs equals the difference between 
the primary beam charge density and Goldreich-Julian charge density. In computing the source 
function of pairs from ICS photons, we consider a hybrid spectrum of RICS and NRICS photons 
from the primary particle, of Lorentz factor 7 , at a given altitude of its acceleration. In the ICS 
photon spectrum, photons with energies e < 2 q B' are produced primarily by RICS. Since the 
RICS spectrum from a single electron cuts off sharply for photon energies e ~ 2 y B', photons 
above this energy, up to a maximum energy of e = 7 will be produced by NRICS. We therefore 
describe the ICS photon spectrum for e < 27 B' by a RICS spectrum based on an expression given 
by Dermer (1990). Assuming delta-function distributions of electrons of Lorentz factor 7 , and 
thermal photons of energy eo = 2.7 kT/mc 2 and incident angles p- < p < p+ in the lab frame, 
Dermer’s (1990) equation (10) for the distribution of scattered photons with energy e s , integrated 
over scattered angles, — 1 < p s < 1 , becomes 


N R (e s ) 

dtde s 


ca T e B 1 n 0 
2 (P+ ~ P -) /3 2 7 3 eg 


[Jo + Ji + -(J2 — Ji)], 


( 20 ) 


where J n = I n (p +) - I n {p~), 


u + = 7e j B 1 max[e s /27 2 ,eo(l - /3p+\, 
u- = yejj 1 min[ 2 e s , e 0 (l - (3p-)\, 


(21) 


and 


Io(u) 


3 


-eg e 2 'ye 0 e s e B {eo + e s ) 

3 u 3 u 2 


( 4/3765 ) 4 
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£ s( e o + e s + 47 2 e 0 e s ) 


u 


- 2'ye\ (e 0 + e s ) In u + e^ri 


/i(u) = u — (it + 1) 1 — 21n(it+l), 
h(u) = u + ln[(u — l) 2 + a 2 ] + (a -1 + a) tan -1 


(u- 1 ) 


( 22 ) 

(23) 

(24) 


where cry is the Thompson cross section, (3 is the electron velocity in units of the velocity of 
light, 6b = B/B cr is the cyclotron energy in units of me 2 , a = 2afes/3 , a/ is the fine-structure 
constant and no = 20cm _3 T 3 is the density of thermal photons of energy eo- Each scattered 
photon is assumed to be emitted along the direction of the particle momentum (i.e. along the 
local magnetic field). Although the distribution of scattered photon angles depends on energy 
(see Fig. 3 of Dermer 1990), the photons at pair-producting energies scattered by relativistic 
electrons are mostly emitted along the field direction. Since we will be treating scattering of both 
upward-moving electrons and downward-moving positrons, we have 


for electrons: /z_ = fi c 

M+ = 1 

for positrons: //_ = —! 


(25) 


M+ = 


(26) 


in equation (21), where 


fJL c = cos 9 C = 


v / 'Tr 


(27) 


is the largest scattering angle at height h above the NS surface, and ry = r pc <C R is the radius of 
the hot PC. 


The ICS photon spectrum for e > 2y B' is produced by NRICS. We use a spectrum for 
relativistic (but non-magnetic) scattering which is based on an expression of Jones (1968), derived 
originally to apply to scattering of an isotropic distribution of photons by a relativistic electron, 
which we have modified to apply to a semi-isotropic distribution of thermal photons 


N NR (e s ) 

dtde s 


no7r?’QcA± 

7 2 eo 


2qlnq + (1 + 2q)(l - q) + 


(grp ) 2 (i-g) 

(1 + (/To) 2 


(28) 


where 

q = r f P Es = r ° = 2 7eoA±- (29) 

To(l -E s ) 7 

Here 77 is the classical electron radius. For upward-moving electrons, A_ = (1 — /3// c ) and for 
downward-moving positrons, A+ = 2. 

Our method for calculating the pair source function from ICS photons is similar to that 
used in Paper I for CR photons. However, unlike in the case of CR, higher generations of pairs 
















-11 - 


from synchrotron radiation of first-generation pairs are important in the screening, since the 
attenuation length for the sychrotron photons is less than the ICS screening scale. As in Paper I, 
the contribution to the pair source function from first generation pairs is computed by dividing 
the ICS spectrum radiated by the primary particle at each step along its path into discrete energy 
intervals. A representative photon from each ICS energy bin is propagated through the local 
field to determine whether it escapes or produces a pair, in which case the location of the pair is 
recorded (for details of such a calculation, see Paper I and Harding et al. 1997). The pairs for 
each ICS spectral energy interval of width Ae s are then weighted in the source function by the 
“number” of RICS or NRICS photons, n R ’ NR , represented by the test photon in that energy bin, 
estimated by integrating the photon spectrum over the bin 


n R ’ NR (e< 


As N R ’ NR (e' s ) , 

V J es dtde' s € ° 


(30) 


where As is the size of the particle path length. 


The contribution of higher generations of pairs to the source spectrum is computed by 
simulating a pair cascade in the method described by Baring & Harding (2001), except that we do 
not include the possibility of photon splitting in the present calculation. By limiting the surface 
magnetic field we consider to Bq < B cr , we can safely neglect splitting. As described in Baring 
& Harding (2001), the created pairs from the first generation make transitions between Landau 
states to radiate synchrotron/cyclotron photons until they reach the ground Landau state. Each 
synchrotron/cyclotron photon is individually traced through the local magnetic field until it either 
escapes or creates a pair, which radiates another generation of photons. A recursive routine is 
called upon the radiation of each photon, so that a large number of pair generations may be 
simulated. 


One notable difference was introduced in the present cascade simulation in order to treat the 
synchrotron radiation from ms pulsars. The low surface fields of the ms pulsars, combined with 
the photon energies e ~ 10 6 — 10' required to produce pairs, give very large Landau states for the 
created pairs, making it impossible to treat the synchrotron/cyclotron photon emission discretely. 
For local fields B < 0.01H cr , we treated the synchrotron radiation from the pairs by the same 
method as for the ICS radiation from the primaries, i.e. as a continuous spectrum. To describe the 
total spectrum radiated by a particle as it decays from a high initial Landau state (i.e. large initial 
pitch angle) to the ground state (zero pitch angle) we use the calculated synchrotron spectrum of 
Tademaru (1973) 

N s (e)de = i(H'sin^ 0 )" 1/2 e- 3/2 de, (31) 

where ip o is the initial pitch angle of the particle and B' = B/B cr is the local magnetic field 
strength. Surprisingly, for the ms pulsars we are required to limit the upper end of the synchrotron 
spectrum at an energy equal to the particle Lorentz factor in order to avoid violating energy 
conservation. This is a well-known problem in high magnetic fields (e.g. Harding & Preece 1987), 
but in this case the very high Lorentz factors of the pairs give a synchrotron critical energy 
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e cr = (3/2)7 2 iT sin'i/’o > 7. Imposing the energy conservation condition limits the number of pair 
generations. The discrete QED treatment of the synchrotron emission in high magnetic fields 
(described above) is able to treat the effect of pair production in the ground Landau state, which 
also severely limits the number of pair generations in fields B ~ 0.2 B cr (Baring & Harding 2001). 


3.2.1. Screening effect at the upper pair front 

Since the ICS photon spectrum has a much weaker dependence on the particle Lorentz factor 
than the CR photon spectrum, the ICS pair source function grows much more slowly with height 
above the PFF. In the case of NRICS, the photon emission rate decreases as ~ I/7, while the 
RICS emission rate first increases sharply at low 7 as the soft photons come into resonance in the 
particle rest frame, and then decreases as I/7 2 , whereas in the case of CR the photon emission 
rate increases as 7. Even though the characteristic ICS photon energies are increasing with 7, the 
stronger decrease in photon number with height as the particle accelerates, combined with the 
decrease in soft photon density and angle with altitude, makes screening by pairs produced by ICS 
photons much more difficult and in some cases impossible. In nearly all cases, the screening scale 
lengths are comparable to the particle acceleration length, requiring different treatment of the 
self-consistent screening calculation from what was used in Paper I for CR screening. The onset 
of pair production and possible screening close to the NS surface couples the entire £) solution 
to the screening scale. We therefore must correct the accelerating Fy as well as the Eu in the 
screening region above the PFF on every iteration of the screening scale length. This is done by 
adjusting the upper boundary of the solution (eqs. [A7] - [A 10] of Paper I) to zo + A s , so that 
the electric potential distribution below the pair front is accurate. We still assume the exponential 
form with scale height (eq. [26] of Paper I) for E\\ above the PFF. But unlike our calculation 
of CR pair screening, the location of the PFF will change and adjust on every interation of A s . 
The exponentional function E^{z ~ zq) = E^ cc (zo) exp[— {z — zo)/A s ] describing the E\\ above 
the PFF is also renormalized at the PFF at each iteration so that the value of E?S c (zq) is equal 
to the maximum value of the E\\ solution between z = 0 and z = zo + A s . We also must take 
into account the increase in the charge deficit in the screening region, which will also cause ICS 
screening to be more difficult. 

Figure 2 shows an example of an ICS pair source function, integrated over energy, so that 
the growth of pairs as a function of distance x = z — zq above the PFF is displayed. The first 
generation pairs (those produced by primary electrons) is almost flat before beginning to decrease. 
Including pairs from higher cascade generations produces an initially sharp increase in number of 
pairs vs. height, but the pair number quickly levels off and also eventually decreases. The ICS 
pair source function is a sharp contrast to the CR pair source function shown in Figure 1 of Paper 
I, which increases nearly exponentially with height above the PFF. This raises the interesting 
question of whether ICS pairs are capable of screening E^ to all altitudes, even if they are able to 
initially screen the local Eu. Figure 3 shows an example of a self-consistent solution for the charge 
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density p(x)/pcj due to polarization of pairs and A s in the screening region. The charge density 
initially increases faster than the charge deficit, which is increasing approximately as x. In this 
case, screening of Ty is achieved locally, until the pair source function begins to decline with x, 
causing the charge density due to pairs to decrease. As the charge density due to pairs decreases 
while the charge deficit is increasing, screening cannot be maintained and the En will again begin 
to grow. We cannot model this effect in the present calculation, which would require the solution 
of Poisson’s equation throughout the screening region. However, it is clear that ICS screening 
cannot short-out En and produce a limit to particle accleration at low altitudes. ICS screening, 
even if locally complete, will only slow the particle acceleration until the ICS pair density falls 
off and beyond this point the particles will resume their acceleration until they can produce CR 
pairs, which will then screen the En at higher altitudes. 

For many pulsars, however, ICS cannot even achieve locally complete screening. Figure 
4 shows computed boundaries of ICS local screening for different assumed PC temperatures, 

Tg = T/10 6 K, as a function of pulsar surface magnetic field and period. Pulsars above and to 
the left of these boundaries can achieve local screening of En at a colatitude of £ = 0.7, while 
pulsars below and to the right of the boundaries cannot. Although the PFF boundaries were 
nearly independent of Tg, the screening boundaries are very sensitive to the temperature of the 
soft photons. For low Tg, only very high-field pulsars, those whose pairs source functions are 
dominated by RICS, can screen with ICS photons. As Tg increases, the boundaries move down, 
but are limited in their movement to the right as they approach the ICS PFF boundaries (Figure 
1). There is also some dependence of the screening boundaries on £, in that the boundary for a 
given Tg will be higher at a lower value of £. This means that screening is more difficult in the 
inner parts of the PC, so that some pulsars may have local ICS screening only in the outer part of 
their PC. 


Even though complete ICS screening is not achieved for all pulsars, there will still be positrons 
which turn around and return to heat the PC. In fact this heating can be substantial since in the 
absence of complete screening which shuts of the acceleration, most of the positrons from ICS 
pairs will return to the PC. Figure 5 shows self-consistent calculations of the returning positron 
fraction, as well as the PFF altitude zo which we have discussed in Section (2.2), as a function 
of £. Figures 5a and 5b, which have the same P and B values but temperatures Tg = 0.5 and 
Tg = 3.0 respectively, illustrate that although there is very little difference in the PFF altitude 
with a large increase in PC temperature there is a large increase in returning positron fraction. 
The dependence of p + /pcj on Tg is not reflected in the analytic estimates (eqs. [18], |19|) which 
are based only on a constant multiple of the charge deficit at the PFF altitude, proportional only 
to zo- The size of p + /p gj is actually dependent on the number of pairs produced above zo which 
is not dependent on the charge deficit if screening is not achieved. Comparison of Figures 5a 
and 5c shows that pulsars with larger periods have relatively larger zo and also larger fractions 
of returning positrons. This is due to the fact that longer period pulsars have more difficulty 
producing pairs due to the smaller PC size and therefore larger radii of curvature of the last open 
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field line. As illustrated in Figure 5d, ms pulsars have very large zq and large fractions of returning 
positrons, but smaller than expected from the analytic solutions, because they do not produce 
enough pairs to screen the field. 

Figure 6 shows the dependence of returning positron fraction on PC temperature Tq and 
surface field strength. The steep dependence of p + / pqj on Tq at low temperatures Tq < 1.0 is due 
to the increase in numbers of high energy ICS photons that produce pairs. For Tq > 1.0, in the 
case of normal pulsars there is local screening at these field strengths, which limits the number of 
returning positrons to be proportional to (zq + A s ) rather than to Tq. The screening therefore 
produces a saturation in the PC heating, which will be seen in the calculation of PC heating 
luminosities to be presented in Section V. 


3.2.2. Screening at lower pair fronts 

As discussed in HM98, the positrons returning from the upper PFF will be accelerated 
through the same voltage drop as the primary electrons and thus will reach the same high energies 
as the primaries that produced them. They will therefore initiate pair cascades as they approach 
the NS surface. HM98 speculated that the pairs from these downward moving cascades may screen 
the E || above the surface since they found that the PFFs of the positrons are located at significant 
altitudes above the surface, due to asymmetry in the ICS between upward moving and downward 
moving particles. Now that we can calculate the fraction of positrons returning from the upper 
ICS PFF, we can investigate the question of screening at the lower PFF. We have simulated 
cascades below the PFFs of test positrons returning from the upper PFF. The source functions 
of pairs from the downward cascades and the charge density are computed in the same way as 
at the upper PFF, and the screening scale height is determined self-consistently. But there are 
several critical differences between screening at upper and lower PFFs. First, the flux of returning 
positrons from the ICS PFF is small compared to the primary flux and the charge density must 
be weighted with the value p + /pcj. Second, the value of E» is very small near the NS surface, 
making it harder to turn around the electrons from the pairs. Third, the charge deficit is much 
smaller than at the upper PFF and decreasing toward the stellar surface, being proportional to 
the altitude of the returning positron pair front. This guarantees that there will be a point above 
the NS surface where p(x)/poj will equal the charge deficit. 

We have computed pair source functions from cascades of downward-moving positrons and the 
resulting charge densities below the lower PFF for a range of pulsar parameters. We find that for 
nearly all pulsars, the charge densities do not reach values high enough to screen E y significantly 
above the NS surface. For pulsars having surface fields B<0.lB cr and periods P < 0.5, pair 
multiplicities from the near-surface cascades are not high enough to balance the low fraction of 
returning positrons, primarily because most pairs are produced in the ground Landau state in the 
high near-surface fields. For pulsars having lower surface fields and longer periods, screening can 
occur for some cases just above the surface at the highest PC temperatures. These pulsars are 
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near the ICS screening boundary (Figure 4). Therefore, according to our calculations it seems 
that ICS PFFs will be stable except possibly near the screening boundaries. 


4. Energetics of primary beam and pulsar death line 
4.1. Luminosity of primary beam 

Before we discuss the pulsar death line, it is instructive to calculate the luminosity of the 
primary electron beam at the PFF. 

The luminosity of the primary beam as a function of 2 (<S 1) can be written as (see eq. [73] 
of MH97, and cf. eqs. [61] of Paper I) 


iprim = ac p(z, Q$(z, £)dS, 

Js(z) 


(32) 


where a = ^ 1 — r g /R ( r g is the gravitational radius of a NS of radius R), p is the charge density 
of the electron beam, <f> is the electrostatic potential, dS = r pc dQ^, and dQ^ = ^d^dcj) is an element 
of solid angle in the PC region (see also Paper I for details). Here the integration is over the area 
of a sphere cut by the polar flux tube at altitude 2 . 

Similar to our derivation of formulae (63)-(65) of Paper I, we can write for the luminosity of 
the primary beam 

Pprirn — ,/pri rn Ryot, ■ (33) 

where E ro t is the pulsar spin-down luminosity (= H 4 RqR 6 /6c 3 / 2 (1), where Bo/f(l) is the surface 
value of the magnetic field strength corrected for the general-relativistic red shift; see eq. [66] 
of Paper I for details), and / pr im is the efficiency of converting spin-down luminosity into the 
luminosity of the primary beam. Thus, by comparing the above two equations, we can get the 
following expression for / pr i m 


/prim = 6 ■ 10 _3 P' 1/2 C 


c 

0.6 


if p < p (CR,R,NR) 

if pit UCR,R,NR) 


(34) 


where C is the altitude scaled by the PC radius, r pc . In the derivation of formula 
k = 0.15 and cos y ~ 1. 


we assumed 


Now we should evaluate expression (|34|) at the altitude beyond the PFF where the electric field 
is screened, or at C* = Co + A s R/r pc . For the CR case C* ~ Co, while for the ICS pair fronts (for 
both R and NR regimes; see also eq. |l7|) 


c* 


1.5Co if P<P* 
2Co if P ~ P< 


(. R,NR) 


(35) 
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After substituting £q we get the following explicit expressions for (* 
Curvature radiation 


ACR)_j 30 P 9/7 /P 4/7 = 2 B% 7 t% /U if P<pS CR) , 
\ 430 P 9/4 /Pi 2 = 4 B% 4 Tq /8 


if P ~ P ( f K) . 


(36) 


Resonant ICS 


d R) = 


10 P 7p /B 12 = b\/ 2 G tI / V2 if P < Pi R) , 
30 P 7 / A /B 3p = b]£tI /% if P ~ P* (i7) , 


(37) 


Non-resonant ICS 


Anr)_I 3 P 7 / 6 /B% 3 = 0.3 bI P t 7/U if P<Pi NR \ 
\ 6 P 7/4 /B 12 = 0.2 B% A tI /8 


if p~p! nr \ 


where Tg = r/10 6 yr = 65(P/Pi 2 ) 2 ■ Inserting expressions (31) - (38) into formulae 
get 

Curvature radiation 


(38) 


, (|34]) we 


LS = 10 32 erg/s \ U -° " q/4 ^ 12 = " U % 2 7 

p 1 0.1 P~ 9 / 4 B l2 = 10 B 12 5/ t 6 9/8 


0.3 P~ 27 / 14 B% 7 = 20 S 12 15/14 r 6 27/28 if P < pj: CR) , 


or 


if p~ p* (CR) , 


4Sm = 1Ql6 ( er g/ S ) V2 ^rot 2 


P 1 / 1A B~ 2 /7 if P<Pl CR \ 


/{(JR) = q 01 
-'prim 


,13/14 29/28 
*12 r 6 
5 r3/4 7/8 

O £> 12 Tg 


0.3 P“ 1/4 if P~pi CR \ 
(CR) 


if P~Pi CR \ 


(39) 

(40) 

(41) 


Resonant ICS 


or 


r (R) = in3 i , / 0.3 P- 13 /e = 30 B 12 13/e r 6 13/12 if P<P* (fl) , 

prim * \ 0.05 P- 11/4 S^ 2 = 20 Bf 2 9/4 r 6 _11/8 if P ~ Pj R \ 


t( R ) — ip 15 Iprp'/s'l 1 / 2 F 1/2 f ^ B 12 11 ^ ~‘ r * 

prim S/ rot j 0.2 P~ 3 / 4 B~ 1/2 if P~ 


if P<P* (R) , 
3 12 if P- P* ( 

f 5 _1 / 6 r 11 / 12 if p<p( R ) 

f( R ) _ ln -2 ) - d 12 r 6 > 

prim 1 0.8 p- 1/4 r fi 5/8 if P ~ P* (R) , 


(42) 

(43) 

(44) 
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Non-resonant ICS 


or 


ANR) 

prim 


10 30 erg/s 


0.3 £^ 3 P- 13 / 6 = 30 bJ /2 t 6 13/12 if P < Pi NR) , 
0.1 B V2 P~ 11/4 = 30 B- 2 7/ \~ UlS if P ~ Pi NR) , 


ANR) 

■^prim 


10 15 (erg/s) 1 / 2 Eg 


0.1 B 12 1/3 P~B 6 if P<Pl NR \ 
0.05 P” 3 / 4 if P~Pi NR \ 


ANR) 

•'prim 


f p 1/2 11/12 ifp<p( NR ) 
,-3 ) n l'2 '6 11 ^ ~ ^ > 

| 2BJV? 8 if P ~ P* NR ^. 


(45) 

(46) 

(47) 


The above formulae are just different representations of the luminosity that a primary (electron) 
beam achieves by the moment the E\\ is screened, for each of the radiation processes we discuss. 
Note that the eight known EGRET 7-ray pulsars fit an empirical relation P 7 oc P r( j t . This is 
reproduced by the above expression for the luminosity of the primary beam at the CR PFF if the 
primary electrons are nearly 100% efficient at conversion of their energy to high-energy 7-rays (see 
Zhang & Harding 2000), but not by the luminosity at the ICS pair fronts. But we have argued 
that the primary particle acceleration is not stopped at the ICS PFF, so that the potential drop at 
the CR PFF determines the final acceleration energy and thus the 7-ray luminosity of the pulsar. 


4.2. Derivation of theoretical pulsar death line 

4-2.1. Analytic approach 

In pulsar physics the term “death line” was originally introduced for radio pulsars and 
refers to a line separating the domain (normally in the P—P diagram) favouring pair formation 
from the domain where pair formation does not occur. Most PC models for radio pulsars imply 
therefore that radio emission turns off (see Sturrock 1971; Ruderman & Sutherland 1975; and 
Chen & Ruderman 1993) if the potential drop required for pair production exceeds the maximum 
which can be achieved over the PC. This concept of death line implicitly assumes that the 
pair-formation is a necessary condition for pulsar radio emission, and that pulsars become radio 
quiet after crossing the death lines during their evolution. Some authors argue that this condition 
may not be sufficient, though. For example, HA01 define a pulsar’s death as the condition 
when pair-production is too weak to generate the pair cascade multiplicity required to screen 
the accelerating electric field. We believe that pair-formation is vital for pulsar operation, but 
we suggest that the pairs required for radio emission may not necessarily screen or shut off 
the accelerating electric field (see also Paper I). In addition, we adopt a canonical definition 
of a pulsar’s death as an essentially pairless regime of operation. In general, to write a formal 
criterion specifying the death line one needs to know the characteristic voltage drop in the pulsar 
acceleration region. For the acceleration model we employ in our study this voltage drop can be 
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calculated as 


A4> acc = 10 10 V B 12 p- 5/2 C 


3 C if P<P* 


{CR,R,NR ) 


if p > p (CR,R,NR)' 


(48) 


The corresponding calculated Lorentz-factor of a primary electron accelerating from the PC is 

j(CR,R,NR) 


'Jacc 


mc- 


r A$ acc = 10 4 J B 12 p- 5 / 2 C 



if p > p(CR,R,NR)' 


(49) 


Thus, our criterion for the pulsar death line separating radio-active from radio-quiet pulsars 
translates into 

'"face A 7min- (50) 

However, the above expression for 7 acc , for fixed values of B and P, still depends on a dimensionless 
acceleration altitude, (, and can not be used therefore in criterion ([5(]). It is important to note, 
that in our numerical calculations, at each altitude, we check the instant value of a particle’s 
Lorentz-factor against the above criterion. If it meets this criterion, we pick the corresponding 
bulk pulsar parameters such as B and P and store them as those of the death line. Generally, to 
make our analytic derivation of the death line self-consistent, we need an additional independent 
equation that relates £ with e.g. B and P. Obviously, the formal substitution of values of £* 
calculated in § 2.2 into 7 acc in criterion (|50|) would lead to a trivial identity, simply because in our 
derivation of £* we have already used the corresponding values of 7 m i n . Instead, we suggest using 
equations (34) and (|49|) to eliminate £ to get 


7 acc ~ 6 • 10 7 / prim Hi 2 P 


-2 


(51) 


The advantage of this formula for 7 acc is that it does not discriminate between unsaturated 
and saturated regimes of acceleration of primary electrons and is factorized by the efficiency 
of the pulsar accelerator, / pr i m , which may be generally regarded as a free pulsar parameter, 
characterizing the efficiency of primary acceleration above the pulsar PC. 


Now, let us use formula (plD in criterion (||b|) to get explicit expressions describing the death 
lines for different underlying radiation mechanisms for pair-producing photons. Note that for this 


purpose in criterion (50) we should evaluate 7 acc at / pr im = .fprlm > where is the minimum 


pulsar efficiency allowing pair formation. This completes our formal definition of death lines 
and makes it physically sensible: for the fixed values of B and P, to specify the onset of pair 
formation, one needs to know the minimum energetics a primary electron beam should have to 
enable pair formation, i.e. one needs to know . We must note that, within the framework of 
this approach, the analytic death lines derived in previous studies implicitly assume some fixed 
value for /“££, the same for all pulsars and for all relevant pair-formation mechanisms, and that 
those values of are not necessarily meaningful. On the contrary, our analytic death lines are 
explicitly determined by a pulsar bulk parameter . The resultant analytic death lines (or 
rather parameter spaces with allowed pair formation) in the P—P diagram read 
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Curvature radiation 


Resonant ICS 


log P > 


f log P - 1-75 log /£& - 16 if P < P ^ R) , 
| log P - 2 log - 17 if P ~ P* (CR) 


(52) 


logP > 


Non-resonant ICS 


| log P - log - 17.4 if P < P* (R) , 
| log P - 1.33 log - 18.9 if P Z Pi R) 


(53) 


logP > 


(NR) 

log P - 2 log - 21.6 if P- Pi 


log P - 1.5 log /“P - 19.8 


if P < P* 


(54) 


/ prim 


The fitting of observational data with the above expressions requires specification of the 
pair-producing efficiency (/“PJ of a pulsar PC accelerator. In this paper we calculate pulsar 
death lines numerically, and by comparing the analytic death lines with our numerical death lines 
we can estimate the values of (see next Section). It is of fundamental importance to know 
if this value of (i.e. / pr i m calculated along a pair death line) is the same for all pulsars, and 
how it depends on the radiation process responsible for the pair formation. We find that the 
analytic death lines defined by (52)-(^4j) with the appropriately chosen value of different 

for each radiation mechanism, fit satisfactory the observational data and are in a good agreement 
with our numerically calculated death lines. This means that our introduction of the parameter 
/prim an d the assumption that this parameter only weakly depends on pulsar B and P values, for 
a given mechanism of generation of pair-producing photons, may be quite justified for the analytic 
tackling of pulsar death lines. 


The analytic expressions above for the pair death lines differ significantly from those derived 
by Zhang et al. (2000). The pair death lines of Zhang et al. (2000) were derived from the condition 
(in our notation) that A4> acc = A<3? max , which implicitly assumes that = 1, whereas we 

allow to be a free parameter. As we shall see in the next section, comparison with numerical 
calculations implies <C 1. Also, Zhang et al. (2000) computed the altitude of the PFF 

and thus the value of A<E> acc from the condition 5 p (e m i n ) = 5 e (7mi n ), that photon attenuation 
length is equal to the length required for the primary electron to produce one ICS photon, since 
the acceleration length S acc <C S p . This introduced a dependence of the PFF altitude on PC 
temperature, which they took to be the self-consistent temperature from PC heating. By contrast, 
our analytic expression for the PFF altitude (eq. [0]) does not depend on PC temperature. Finally, 
Zhang et al. (2000) included pairs from only RICS and did not include pairs from NRICS in their 
ICS death lines. 
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4-2.2. Comparison with numerical calculations 


The discussion and results of our numerical calculation of the CR and ICS pair death lines 
were presented in Figure 1. Since our numerical solution for the location of the PFFs are computed 
by iteration, we can unambiguously determine the pair death lines without having to define the 
efficiency, , as was needed in the analytic expressions above. 

We must note that in our derivation of the above analytic death-line conditions we did not 
require the value of accelerating potential drop to be maximum at the PFF. In other words, we did 
not impose any ad hoc constraint on the acceleration altitude, thus making the latter consistent 
with the accelerating potential drop. It is remarkable that for each of the radiation mechanisms 
the parameter (minimum efficiency of acceleration of primary beam required to set pair 

formation) remains nearly constant along the numerical death line. By comparing our numerically 
computed death lines with those given by (|52|)-(|54|) we find that 7™^ ~ 0.03 — 0.1 for CR, 0.003 - 
0.02 for RICS and 4 — 6 x 10 ~ 3 for NRICS pair fronts. These efficiencies are still small compared 
to the maximum efficiency we esimate in the last paragraph of this Section. Another important 
finding is that production of pairs in all observed pulsars by CR photons would require 1 

indicating the difficulty in interpreting the observational data in terms of a CR mechanism alone. 
On the contrary, ICS-based death lines imply much less consumption of pulsar spin-down power 
for pair creation. This result and our finding that the ICS-generated pairs tend to only partially 
screen the accelerating electric field may suggest the occurrence of “nested” pair formation regions: 
the lower-altitude pairs produced by the ICS photons and the higher-altitude pairs produced by 
the CR photons generated by particles accelerating through the region of the ICS pair formation. 

This fact implies that if a pulsar is below the CR death line in Figure 1, and if any ICS 
screening is inefficient, then the primary beam acceleration is not limited by pair production and 
will be producing ample CR high-energy photons. The upper limit for the pulsar 7 - 1 'ay luminosity 
in this case can be estimated by using equation ( |32|) with the integrand evaluated at the maximum 
value of the potential given by formula (13) of HM98 (see also eq. [A4] of HM98): 


L 


7 , max 


= C 


ls(v) 


max dS( V ), 


(55) 


where $ max ~ (1/2)[TIR/ cf (1)\k&o , 4>o = {QR/c)BqR, and dS{rf) = S(rj)dQ^/TT, where S(r /) is 
the spherical cross-sectional area of the polar flux tube at radial distance r (= rjR). Note that in 


formula (55) the integral c/ apdS represents the total electron current flowing from the PC and is 
a constant. Thus, integrating over £ and q i, we arrive at 


L 


7 , max 


1.57tk(1 — n)E V ot ~ 0.6-F ro t. (56) 

Here we used k = 0.15 and cosx ~ 1. We predict, therefore, that pulsars tend to be efficient 7 -ray 
sources if they are to the right of their CR death lines in B-P diagram. Note that expression (56) 
can also be used to set the possible upper limit for n, which is around 0.3. 
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4.3. Characteristic voltage drop at the pair front 


Now we would like to demonstrate one more remarkable property of a pulsar PC accelerator 
Let us take the expression for the electric potential as a function of ( (see eq. |48|]) and evaluate 
it at the screening altitude (at ( = £*)• After substituting the corresponding expressions for 
(see eqs. fl3 


!^]) into equation (48) we get the following self-consistent formulae for the critical 
voltage drop (characteristic voltage drop at the screening altitude) 


Curvature radiation 


A4w(C*) = 1013 V 


2 (re/P) 1 / 14 if P<P* (aR) , 
p-i/A if P~pf CR) . 


Resonant ICS 


A<h acc (C*) = 10 11 V 


3 P ~ 7 /6 t i / 2 if p<pj R \ 

2 P~ 5 / 4 Tg /4 if P~P* {R) , 


Non-resonant ICS 


A$ acc (C*) = 10 11 V 


2 p-VVg 176 if P<Pi NR \ 
p-3/4 if P~P* (Ari?) . 


(57) 


(58) 


(59) 


These formulae show very weak dependence of critical voltage drop on pulsar parameters B and 
P, especially for the CR case, as also shown by HM98. This fact simply indicates that the critical 
voltage drop needed for the ignition of a pair formation is mainly determined by the pair-formation 
microphysics itself. As a result, as we move from the lower left to the upper right corner of 
the P-P diagram, the variation (by orders of magnitude) in the altitude of the PFF effectively 
compensates the corresponding change in voltage drop due to its scaling with pulsar parameters 
B and P, thus maintaining it near its critical value. 


Note that in similar expressions derived by HA01 (see their eqs. [73]-[75]) the characteristic 
voltage drop varies by ~ four orders of magnitude over the whole range of pulsar spin periods 
~ 0.001 - 10 s. 


5. Polar cap heating luminosities and surface temperatures 
5.1. Analytic estimates 

To estimate the efficiency of PC heating by returning positrons, we can use formulae (64), 
(65) of Paper I again evaluated at z*. 
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Thus, setting k = 0.15 and cosx ~ 1, we get for the heating efficiency, /+ = L + /E rot , 
Resonant ICS 


/f 


-5p-l/2 3/2 


10"°P 


2.5 if P<Pi R \ 

1.5 ifP~P* (i7) , 


(60) 


Non-resonant ICS 


j(NR) 


10 - 5 P 1 / 2 r 6 


1.0 if P< P* ( 

0.6 ifP~P* (JVR) . 


(61) 


Here we shall summarize the explicit expressions for the estimated PC luminosities and surface 
temperatures due to the heating by positrons returning from the pair fronts set by the CR 
and ICS photons. Our estimates of the PC temperature are based on a standard formula 
T pc ~ (L+firr^asB) 1 / 4 , where cfsb is the Stefan-Boltzmann constant. This formula implies 
therefore that the PC is heated homogeneously, and that the heated area is confined by the PC 
radius. 


Curvature radiation 


l (CR) = 10 si erg / g p - 1/2 | 


0.4 P -5 / 14 t 6 

1.0 

1/7 if P<P* (Ci?) , 
if P-PJ CR) , 

(62) 

p—5/56_~ 1/28 
r r 6 

2 

if P<Pj CR) , 
if P~ Pi C7?) , 

(63) 


Resonant ICS 



10 28 erg/s 



p—5/2 1/2 / 0.8 if P<Pi K) , 

6 \ 0.5 if P ~ Pi R) , 

8 x 10 3 K P“ 9/4 r 6 1/8 , 


(64) 

(65) 


Non-resonant ICS 


l {NR) = 10 28 erg/s p- 3/2 


0.3 if P<Pl NR \ 
0.2 if P~Pl NR \ 


( 66 ) 


T^c R) ~ 4 x 10 4 K P~ 5/4 . 


(67) 


We caution that the above formulae are upper limits to the PC heating luminosity, in that they 
assume locally complete screening. As we have shown in Section 3, such screening does not occur 
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in many cases, most notably for older pulsars. Thus as will be shown by the numerical calculation 
presented in the next Section, eqs. ( | 66 | ) and (67) do not accurately predict L + and T pc for ms 
pulsars, but give values that are much too high. 


Let us compare the above temperature estimates with those of HA01. At B 12 = 1 and 
P = 0.1 s our estimate of T pc R is roughly the same as the corresponding estimate presented 
by HA01 (see their eqs. [58]). For the ICS cases (see their eqs. [56], [57] and [59], [60]) their 
estimates are systematically bigger, by a factor of 3 to 12. However, for the ms pulsars we predict 
higher temperatures. Also, the effect that the ICS contribution to the PC temperatures dominates 
over the CR contribution for ms pulsars is more pronounced in our formulae. In addition, for 
long-period pulsars our formulae predict more drastic decline in the PC heating temperatures with 
a pulsar period. We also comment that both the analytic estimates and numerical calculations of 
L + depend on the value of n, related to the compactness and moment of inertia of the NS. 


5.2. Numerical calculations 


We numerically evaluate the luminosity of PC heating due to ICS PFFs in a similar way to 
that given by equation (62) of Paper I. The fraction of returning positrons times the potential 
drop between the PFF and the surface is integrated in £ across the PC. However, since in the 
case of ICS PFFs the potential varies significantly across the width of the PFF, we also integrate 
the product of the potential and the returning positron fraction at each altitude between zo and 
zq + A s : 

L + = 2acS pc rjlj^ p + (z - z 0 ,£,)<S>(z,£) dz £d£ = L+(£) £d£ (68) 

where S pc = 7 rHR 3 /c/(l) is the area of the PC. Figure 7 shows the dependence of L + (£) on £, 
which reflects the distribution of heating across the PC. Most of the PC heating due to the ICS 
pair front occurs at small £, near the magnetic pole for normal pulsars, but in the outer part of 
the PC for ms pulsars. 


The total positron heating luminosity, scaled with the spin-down luminosity, as a function 
of characteristic pulsar age, r = P/2P is shown in Figure 8 for two different values of PC 
temperature. The numerically computed L + /E mt increases with r, in agreement with the analytic 
estimate (eqs. 


and [61]), as long as there is locally complete screening. When complete 
screening is no longer achieved, at a t value that depends on P and Tg, the heating rate saturates 
because the number of pairs produced is not sufficient to return a positron flux proportional to zq. 
Even at very high PC temperatures, the ms pulsars have reached saturation of the PC heating 
rate because they are beyond the ICS screening boundary. The heating luminosity is larger for a 
lower value of Tq before saturation because the PFF is at higher altitude and the potential drop 
is higher, giving the positrons more energy before they reach the NS surface. Figure 9 displays 
the same calculations of the PC heating luminosity, not scaled to the spin-down luminosity, to 
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more easily compare with observations. In these Figures we have also plotted for comparison the 
calculations of PC heating luminosity from the CR PFFs presented in Paper I. The high r end of 
the CR heating rate curves marks the CR PFF pair death line for that period. Where heating from 
positrons returning from the CR PFF is present, it is several orders of magnitude larger than the 
heating from the ICS PFF, so that in these pulsars heating from the ICS pairs makes a negligible 
contribution to the total PC heating rate. For normal pulsars, heating from ICS pair fronts is 
not high enough to be detectable at present, even in the absence of CR pair heating. However, 
emission from ICS pairs provides the only source of PC heating for the known ms pulsars, which 
have ages r > 10 s yr, since they cannot produce CR pairs. This emission may be detectable if the 
PC temperatures are above 10 6 K. The relatively high ICS pair heating rates in the ms pulsars is 
due primarily to the higher voltage drop necessary to make pairs. Even though the fraction of 
returning positrons is not large, the positrons gain more energy before hitting the NS surface. 


The analytic estimates for L + are in reasonably good agreement with our numerical 
calculations (within a factor of 10 ) in the regime where locally complete screening is occurring. 
Where compete screening is not achieved, the analytic estimates of L + and also of T pc given in 
Section 5.1 will greatly overestimate the true values. This is especially true in the case of the ms 


pulsars. We caution that equations (54) 
screening boundaries of Figure 4. 


(67) cannot be used for pulsars that are beyond the 


In Figure 10, we show the predicted luminosities from ICS PC heating, L + as a function 
of PC surface temperature (solid lines) for two different pulsar periods and a surface field of 
Bq = 4.4 x 10 8 G, so that the corresponding characteristic ages would be r = 1.4 x 10 9 yr for P 
= 2 ms and r = 8.75 x 10 9 yr for P = 5 ms. There is a definite dependence L + oc T 2 seen in 
these curves, which cannot be modeled analytically since these cases do not have full screening. 
Also plotted (dashed lines) are the luminosities, Lbb = emitted by a blackbody at PC 

radiating temperature T and different heated surface areas A. The intersections of the curves 
roughly indicate values of temperature where self-sustained heating, (i.e. where the surface 
emission at a given temperature supplies just enough returning positrons from ICS pairs to 
maintain the PC at that temperature) is possible for heated PCs of a particular area. Since the 
standard PC area is A pc = 3.3 x 10 11 cm 2 for P = 2 ms and A pc = 1.3 x 10 11 cm 2 for P = 5 ms, 
self-sustained PC heating emission requires heated and radiating areas much smaller than the 
entire PC area. The intersection points of the blackbody and L + curves in Figure 10 are not 
entirely self-consistent self-sustained heating models since we assume a radiating area of A pc for 
our calculations. The fully self-consistent models would have radiating areas somewhat larger at a 
given temperature, but still smaller than A pc . 

Of the six millisecond radio pulsars which have been detected as pulsed X-ray sources, 
most have narrow pulses and power-law spectra indicating that their emission is dominated 
by non-thermal radiation processes. However, two pulsars, PSR J0437-4715 (P = 5.75 ms, 
r = 4.6 x 10 9 yr) and PSR J2124-3358 (P = 4.93 ms, r = 7.3 x 10 9 yr), have possible thermal 
emission components which would imply that some surface heating is taking place. The emission 
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from PSR J0437-4715 has a dominant power-law component but a two-component model is 
required for an acceptable fit. A two-component power-law plus blackbody fit to combined EUVE 
and ROSAT data (Halpern et al. 1996) give a temperature of T = (1.0 — 3.3) x 10 6 K, luminosity 
L = 8.4 x 10 29 ergs ” 1 and emitting area A = 7.8 x 10' — 1.1 x 10 10 cm 2 for the thermal component. 
Recent Chandra observations of PSR J0437-4715 (Zavlin et al. 2001) confirm that at least one 
thermal component plus a power law is needed to fit the spectrum. Their preferred model consists 
of a two-temperature thermal blackbody, with a hotter PC core of T core = 2.1 x 10 6 K and 
Rcore — 0-12 km and a cooler PC rim of T r j m = 5.4 x 10 5 K and R r im = 2.0 km. A blackbody fit 
to ASCA emission from PSR J2124-3358 (Sakurai et al. 2001) yields a temperature 3.6;tg' 9 g x 10 6 
K, luminosity L = 4.8 x 10 29 ergs _1 and PC emitting area A = 1.4lg'g x 10 7 cm 2 . 

Comparing our results in Figure 10 with the observed values of T, L and A for PSR 
J0437-4715 and PSR J2124-3358, we see that the predicted luminosities, areas and temperatures 
for self-sustained PC heating are in a range comparable to observed values. Thus ICS pair fronts 
could be a plausible source of PC heating for some ms pulsars. 

Recent Chandra high-resolution X-ray observations of the galactic globular cluster 47 Tuc 
(Grindlay et al. 2001a) have detected all 15 known ms radio pulsars in the cluster as well as a 
number of other suspected ms pulsars. It is thought that the X-ray source population of 47 Tuc is 
dominated by ms pulsars having soft spectra and luminosities Lx ~ lO'^ergs” 1 . Grindlay et al. 
(2001b) have found a dependence log(L.Y) = —(0.32 ± 0.1) log(r) + 33.3 for the ms pulsars in the 
cluster, not too different from what we have found for PC heating from ICS pair fronts (cf. Figure 
9b). 


6. SUMMARY AND CONCLUSIONS 

We have explored production of electron-positron pairs by photons produced through ICS 
of thermal X-rays by accelerated electrons above a pulsar PC. Since the accelerating primary 
electrons can produce pairs from ICS photons at much lower energies than are required to produce 
pairs from CR photons, it is very important to investigate the consequences of ICS pair fronts 
for A| screening and for PC heating. We have defined “pair death” lines in Bq — P space as 
the boundary of pair production for pulsars, having dipole magnetic fields of surface strength 
6.4 x 10 19 (RP ) 1 / 2 G. Operationally, the existence of a pair front is determined by a finite solution 
to equation (0) for Sq, the altitude of the onset of pair creation. Although we are able to give 
analytic formulae for the altitude of the PFFs for different radiation processes, the location of the 
pair death lines must be determined numerically. The existence of a pair front requires much less 
than one pair per primary electron (but still many pairs from the whole PC beam), since the very 
first pairs are created in the declining high-energy tail of the radiated spectrum, and also much 
less than the number of pairs per primary required for screening of Eh. We find that virtually all 
known radio pulsars are capable of producing pair fronts with ICS photons. A smaller number, 
less than half, are able to also produce pairs via CR. If the acceleration model we use is correct 
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and pulsars have dipole fields, then this result implies that relatively few pairs are required for 
coherent radio emission. 

Self-consistent calculations show that ICS pair fronts produce lower fluxes of returning 
positrons and lower PC heating luminosities than CR pair fronts. This is due to the higher 
efficiency of the ICS process in producing pairs at lower altitudes where both the charge deficit 
required to screen (and thus the returning positron flux) and the accelerating voltage drop are 
smaller. For pulsars with surface magnetic fields in the “normal” range of 10 11 — 10 13 G, ICS 
heating luminosities are several orders of magnitude lower than CR heating luminosities. However, 
for ms pulsars having surface fields in the range 10 8 — 10 10 G, production of any pairs requires such 
high photon energies that ICS pair fronts occur at higher altitudes, where acceleration voltage 
drops are high enough to produce significantly more PC heating. Since most ms pulsars cannot 
produce pairs through CR, ICS pair fronts provide the only means of external PC heating. We find 
that for surface temperatures T ~ 10° K, ICS heating luminosities are in the range of detection. 

We find that ICS pairs are able to screen the local En in some pulsars having a high enough PC 
temperature, but that this local screening will not produce a complete screening of the accelerating 
field at all altitudes and thus will not stop the acceleration of the primary beam. This is because 
the number of ICS pairs grows slowly, on a scale length comparable to the altitude of the PFF, 
and then declines while the charge deficit which maintains E y continues to increase with altitude. 
Even if the ICS pair production is vigorous enough to achieve local screening of E\\, it eventually 
cannot produce the charge density (from returning positrons) to keep up with the increasing 
charge deficit (produced by the combination of flaring field lines and inertial frame-dragging). The 
primary particles may slow their acceleration briefly, due to the local screening, but will resume 
acceleration once the ICS pair production declines. At higher altitudes, those pulsars to the left 
of the CR pair death line will reach the Lorentz factors (7 ~ 2 x 10') required to produce CR 
pairs. In contrast to ICS pair fronts, the growth of pairs above the CR PFF is rapid and robust 
due to the sensitivity of CR photon energy and emission rate on particle Lorentz factor, producing 
complete screening of E\\ in a very short distance (cf. Paper I). In pulsars to the right of the CR 
pair death line, there is not complete screening and the primary particles will continue accelerating 
to high altitude with their Lorentz factor being possibly limited by CR reaction. 

We have also investigated the proposal by HM98 that pairs produced as the returning 
positrons are accelerated toward the NS may be able to screen the E\\ above the surface. Using our 
calculated values of the returning positron flux in cases where local screening has been achieved at 
the ICS pair fronts, we find that in most cases screening does not occur at a significant distance 
above the NS surface to cause disruption of a steady state or formation of pair fronts. In the cases 
where screening does occur near PFFs significantly above the surface, the pulsars are near the ICS 
screening boundary. The resulting instability would then not move the start of the acceleration to 
higher altitudes, as HM98 had envisioned, but would probably weaken or disrupt the screening at 
the upper ICS pair front, resulting in a decrease in returning positron flux which would weaken 
the screening at the lower PFF, etc. 
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PC heating by CR pair fronts will dominate for pulsars to the left of the CR pair death 
line, while ICS pair fronts will supply the PC heating for pulsars to the right of the CR pair 
death line. While we have given analytic expressions for the fraction of returning positrons 
and PC heating luminosities from ICS pair fronts, these are only good estimates above the ICS 
screening boundaries of Figure 4. Below the ICS screening boundaries, where local screening is 
not achieved for that PC temperature, the numerical values of returning positron fraction and 
heating luminosity fall well below the analytic estimates. This will be true for ms pulsars, nearly 
all of which are below the ICS screening boundary for a PC temperature of even T = 10 7 K. 

Our results are dependent on a number of assumptions inherent in our calculations. First, 
we have assumed that the pattern of thermal X-ray emission is from a heated PC and is a pure, 
isotropic blackbody. According to recent calculations (Zavlin et al. 1995) of radiation transfer 
in magnetized NS atmospheres, the thermal emission is not pure blackbody or isotropic, but a 
somewhat cooler blackbody consisting of pencil and fan beam components. Both effects of full 
surface (cooler) emission and pencil beaming would tend to decrease ICS radiation and thus ICS 
pair production efficiency. However, a large fan beam component would tend to increase ICS 
efficiency. Second, we have used a hybrid scheme to describe the ICS radiation spectrum in which 
the RICS has been treated as classical magnetic Thompson scattering and NR (or continuum) 
ICS has been treated as relativistic but non-magnetic. In reality, both RICS and NRICS should 
be treated as a single process with one cross section. While the magnetic QED scattering cross 
section has been studied for some time (e.g. Herold 1979, Daugherty & Harding 1986), simple 
expressions in limited cases are only beginning to become available (e.g. Gonthier et al. 2000). 
Our present treatment is probably accurate for magnetic fields B<0.2B C1 - which includes most of 
the radio pulsars. 

The location of a pulsar relative to the pair death lines may be important not only to its radio 
and thermal X-ray emission characteristics, but also to its high-energy emission properties. As we 
have argued in this paper, ICS pair fronts will not limit the acceleration voltage drop in pulsars 
but that acceleration will continue until it is limited by a CR pair front. The voltage drop at the 
CR pair front (together with the size of the PC current) is therefore expected to determine the 
high-energy emission luminosity. In Section 4, we have noted that the acceleration voltage drop of 
pulsars that produce CR pair fronts is remarkably insensitive to pulsar parameters, leading to the 
prediction that high-energy luminosity, Lhe, should be simply proportional to PC current (which 
is proportional to ), which seems to be borne out by observations (e.g. Thompson 2000). 
However, pulsars that do not produce CR pair fronts do not have such a limit on acceleration 
voltage drop and should depart from the Lhe oc E roi 1/2 dependence, and approach a Lhe oc E ro t 
dependence. Indeed, such a departure must occur if they are not to exceed 100% efficiency in 
converting rotational energy loss to high-energy emission. We predict that this change in The 
dependence should occur along the CR pair death line. None of the pulsars which have detected 
7 -ray emission are to the right of the CR pair death line although some, such as Geminga, are 
close. The Large Area Gamma-Ray Telescope (GLAST) will have the sensitivity to detect 7 -ray 
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emission from significant numbers of radio pulsars beyond the CR pair death line and so should 
be able to test this prediction. 


We thank Josh Grindlay, Jon Arons, George Pavlov and Bing Zhang for comments and 
discussion. We also thank the referee for valuable comments. This work was supported by the 
NASA Astrophysics Theory Program. 
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CR, ICS Pair Front 



Fig. 1.- Boundaries defining regions in surface magnetic field .Bo, in units of critical field, 
versus pulsar period, P, where pulsars are capable of producing pairs through curvature radiation 
(CR) or inverse-Compton radiation (ICS). The ICS curves are labeled with different values of PC 
temperature, Tq, in units of 10 6 K. Also shown are radio pulsars from the ATNF Pulsar Catalog. 
Pulsars to the right of each line cannot produce pairs by the corresponding process. 
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Fig. 2.— Example of electron-positron pair source function integrated over energy, as a function of 
distance, x, above the pair formation front produced by a pair cascade initiated by inverse-Compton 
radiation from primary electrons. The vertical axis measures the number of pairs produced in each 
spatial bin, normalized per primary electron. Both curves are for P = 0.1 s, B/B cv = 0.1, Tq = 0.5 
and colatitude £ = 0.7, in units of PC half-angle. 
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Fig. 3.— Self-consistent solution for the charge density, p(x), due to trapped positrons which 
asymptotically approaches the charge deficit, Ap , needed to screen the electric field, Eji, above the 
PFF, modeled as a declining exponential with screening scale height, A S R. x is the distance above 
the pair front. The pulsar parameters are the same as those of Figure 2 except that Tq = 1 here. 
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ICS Screening Boundary 



Fig. 4.— Boundaries in surface magnetic field Bq (in units of critical field) versus pulsar period, P, 
above which locally complete screening of E\\ occurs above the ICS pair front. Curves are labeled 
with values of the PC temperature Tq, in units of 10 6 K. The dashed curve is the CR pair boundary 
(as shown in Figure 1). 
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Fig. 5.— Solutions for the returning (trapped) positron density, p + /P g j> normalized to the 
Goldreich Julian density and pair formation front altitude, zq, in units of NS radius, as a function 
of magnetic colatitude, £, which has been normalized to the PC half-angle. 
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Fig. 6.— Solutions for the returning (trapped) positron density, p + /P gj , normalized to the 
Goldreich Julian density as a function of PC temperature Tq, in units of 10 6 K, for different 
values of surface magnetic field strength in units of the critical field, B/B cr . 
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Fig. 8.— PC heating luminosity from ICS pair fronts, L + , normalized to the spin-down energy 
loss rate, E rot , as a function of the characteristic spin-down age, r = P/2P, for different pulsar 
periods, as labeled, and PC temperature a) Tq = 1.0 b) Tq = 3.0. Closed symbols designate locally 
complete screening and open symbols indicate that no screening occurs above the ICS pair front. 
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Fig. 9.— PC heating luminosity, L + , as a function of the characteristic spin-down age, r = P/2P, 
for different pulsar periods, as labeled, and PC temperature a) Tq = 1.0 b) Tq = 3.0. Solid 
curves show heating luminosities from ICS pair fronts and dashed curves show heating luminosities 
from CR pair fronts (from Paper I). Closed symbols designate locally complete screening and open 
symbols indicate that no screening occurs above the ICS pair front. 
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Fig. 10.- PC heating luminosity L + (solid curves) for different pulsar periods, as labeled, 
and blackbody luminosity Lbb = AusbT 4 for different areas, A as labeled, as a function of PC 
temperature, T. 




